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We present first principle calculations of the structural and electronic properties of the CuS covel- 
lite material. Symmetry-lowering structural transition is well reproduced. However, the microscopic 
origin of the transition is unclear. The calculations firmly establish that the so far controversial Cu 
valency in this compound is 1.33. We also argue that recently reported high-temperature supercon- 
ductivity in CuS is unlikely to occur in the stoichiometric defect-free material, since the determined 
Cu valency is too close to 1 to ensure proximity to a Mott-Hubbard state and superexchange spin 
fluctuations of considerable strength. On the other hand, one can imagine a related system with 
more holes per Cu in the same structural motif {e.g., due to defects or O impurities) in which case 
combination of superexchange and an enlarged compared to CuS Fermi surface may lead to uncon- 
ventional superconductivity, similar to HTSC cuprate, but, unlike them, of an /-wave symmetry. 
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Copper sulfide in the so-called covellite structure (Fig. 
[1]) has recently attracted attention due to a new re- 
port about possible superconductivity at 40 K^. This 
report has been met with understandable skepticism, 
because previous researches'*^ reported reproducible su- 
perconductivity at rather low temperatures, around 1-2 
K. On the other hand, inspection of the literature re- 
veals that reported physical properties of covellite are 
drastically different in different papers. For instance, 
one paper reported a well defined Curie- Weiss magnetic 
susceptibility^, while others observed a nearly constant 
behavior consistent with the Pauli susceptibility in ab- 
sence of any local moments. 




FIG. 1: (color online) Crystal structure of covellite in the 
high-symmetry phase. Large yellow spheres indicate sulfure, 
and blue spheres copper. The dark spheres form the pla- 
nar CuS layers (Cul), and the light spheres form the warped 
CU2S2 bilayers (Cu2). Fat yellow sticks indicate strong cova- 
lent bonds inside the S2-S2 dumbells. 



Structural properties of the covellite are also intrigu- 
ing. At room temperature it consists of triangular layers 
of Cu and S, stacked as follows, using standard hexagonal 
stacking notation: Cul and SI form layers A and B, at 
the same height, so that Cul has coordination of three 
and no direct overlap. Cu2 and S2 form layers B and C, 
so that Cu2 is directly above SI and bonds with it, too, 
albeit more weakly than to S2. Thus, compared to the 
S2 layer, the Cu2 layer is closer to the Cul-f-Sl one, and 
Cu2 appears to be inside a tetrahedron, closer to its base. 
The next layer is again C, so that two S2 atoms are right 
on top of each other and form a strongly covalent bond, 
the shortest bond in this system, essentially making up 
an S2 molecule. 

At T = 55 K the system spotaneously undergoes a 
transition from a hexagonal structure to a lower symme- 
try orthorhombic structure. To a good approximation, 
the transition amounts to sliding the Cu2-S2 plane with 
respect to the Cul-Sl plane by 0.2 A, and the two neigh- 
boring Cu2-S2 planes by 0.1 A with respect to each other, 
in the same direction. The bond lengths change very lit- 
tle, one of the three Cu2-S2 bonds shortens by 0.04 A, 
and the S2-S2 bonds lengthens by 0.05 A, and all other 
bonds remain essentially unchanged. Note that such 
transitions are quite uncommon for metals, but rather 
characteristic of insulating Jahn- Teller systems. Trans- 
port properties are hardly sensitive to this transition, 
which is however clearly seen in the specific heat. 

Thus there are three questions to be asked. First, what 
is the nature of the low temperature symmetry-lowering? 
Second, why some experiments indicate pure Pauli sus- 
ceptibility, while others observe local moments (through 
Curie- Weiss behavior)? Third, why one particular exper- 
iment sees indications of high temperature superconduc- 
tivity, while others do not? Of, course, there is always a 
chance the "outliers" experiments are simply incorrect, 
but it is always worth asking the question, whether some 
sample issues may 

possibly account for such discrepancies. 

In order to adress the first question, we have performed 
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TABLE I: The calculated total energy (meV/cell) of the low- 
temperature orthorhombic and the high-temperature hexago- 
nal structure, using either the experimental or the calculated 
optimized parameters. Structural parameters used, as well 
as selected bond length (A) are also shown. Note that one 
unit cell includes 6 formula units. The cell volume is given in 
The last column corresponds to the orthorhombic struc- 
ture with internal coordinates optimized, while keeping the 
experimental unit cell. 





H-exp 


O-exp 


H-calc 


O-calc 


O-c.o. 
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3.789 


3.760 


3.807 


3.793 


3.760 
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3.789 


6.564 


3.807 


6.623 


6.564 
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16.321 


16.235 


16.496 


16.453 


16.235 


ZCu2 


0.1072 


0.1070 


0.1069 


0.1077 


0.1083 


ZS2 


0.0611 


0.0627 


0.0639 


0.0646 


0.0651 


ycul 


n/a 


0.6377 


n/a 


0.6227 


0.6077 


ycu2 


n/a 


0.3372 


n/a 


0.3410 


0.3413 


ysi 


n/a 


0.3068 


n/a 


0.2917 


0.2760 


ys2 


n/a 


0.0008 


n/a 


0.0064 


0.0069 


Cul-Sl 


3x2.19 


2x2.18 


3x2.20 


2x2.20 


2x2.20 






2.17 




2.19 


2.19 


Cu2-Sl 


2.33 


2.33 


2.36 


2.33 


2.34 


Cu2-S2 


3x2.31 


2x2.30 


3x2.31 


2x2.30 


2x2.30 






2.28 




2.33 


2.33 


S2-S2 


1.99 


2.04 


2.11 


2.13 


2.12 


Volume 


202.9 


200.3 


207.0 


206.7 


207.3 


Energy 





-85 


-258 


-265 


-189 



density functional calculations (DFT) of both hexago- 
nal (H) and orthorhombic (O) structure. First, we op- 
timized the crystal structure using the standard VASP 
program with default settings (including gradient correc- 
tions), and starting from the experimental structure as 
reported in Ref. H. 

After that, all calculations in the determined crystal 
structures were performed using the standard all-electron 
LAPW code WIEN2k. We have also verified that the 
calculated forces in the optimized structures are small 
enough. As a technical note, to obtain full convergences 
in the energy differences we had to go up to -R-ftTmax = 9. 

The results are shown in Table 1 . Even though there is 
some discrepancy between the calculated and the exper- 
imental low-temperature structures (mostly in terms of 
an overall overestimation of the equilibrium volume), the 
correct symmetry lowering is well reproduced. In fact, 
given that only one paper has reported internal positions 
for the orthorhombic structre, and the same paper found 
a Curie- Weiss law, suggesting, as discussed below, crys- 
tallographic defects in their sample, it is fairly possible 
that the calculations predict the structure of an ideal 
material better than this one experiment has measured. 

A more important question now is, what the mech- 
anism for this well-reproduce symmetry lowering can 
be? Ionic symmetry-lowering mechanisms (such as Jahn- 
Teller) are excluded in a wide band metal like CuS. Typ- 



ically, a lower symmetry is stabilized in a metal if it re- 
sults in a reduced density of states at the Fermi level 
( "quasinesting mechanism"). However, the density of 
states at the Fermi level does not change at this tran- 
sition (Fig. [2]), and the states below Fermi actually shift 
slightly upward. Thus, one-electron energy is not the 
reason for the transition. A look at the calculated Fermi 
surfaces (Fig. 0]) shows that while they become more 2D 
in the orthorhombic structure (the in-plane plasma fre- 
quency remains the same, w 4.0 eV, while that out of 
plane, 1.36 eV, drops by 12%), there is no shrinkage in 
their size. 
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FIG. 2: (color online) Calculated density of states in 
the high-temperature ( "hex" ) and low temperature ( "ortho" ) 
structures, using in both cases optimized parameters. 

We cannot say definitively what causes the low- 
temperature symmetry lowering in CuS, but we can say 
confidently that it is not van der Waals interaction as con- 
jectured in Ref. H (for one, it would not be reproduced 
in LGA/GGA calculations with requested accuracy, and, 
also, as dicussed below, the interplanar bonding is cova- 
lent, and not van der Waals), and not a typical metallic 
mechanism driven by Fermi surface changes. One can- 
didate is ionic Coulomb interaction. Indeed the calcu- 
lated Ewald energy is noticeably lower in the orthorhom- 
bic structure. However the Ewald energy is only part 
of total electrostatic energy, so from this fact alone one 
cannot derive definitive conclusions. 

Let us now discuss the electronic structure. Since the 
differences between the two structures are very small, we 
shall limit our discussion by the high-temperature hexag- 
onal structure. The calculated band structure is shown 
in Fig. [3l Note two sets of bands, one at -7 eV and 
the other at 1 eV, of strong S2-p^ character. These are 
bonding and antibonding bands of the S2-S2 dumbells. 
Historically, there has been a heated discussion of the 
Cu valency in this compound, and what is an appro- 
priate ionic model. Both (Cu+i)3(S^2)(-s-i) (ReL i) 
and (Cu+i)3(S2"^)(S-2) (Ref. i) have been discussed, 
assuming monovalent copper. On the other hand, XPS^ 
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FIG. 3: Calculated band dispersions in the hexagonal struc- 
ture. The points F, M and K are in the central plane {kz — 0) 
and A, L and H in the basal plane {kz = tt/c). In the top 
panel, the width of the lines is proportional to the amount of 
the S2-pz character in the corresponding states. 



and NQR^ data indicated Cu valency larger than 1, but 
smaller than 1.5. From our calculations it is immediately 
obvious that SI is divalent, while S2 is monovalent (the 
antibonding band of the S2-S2 dimer is 1 eV above 
the Fermi level, while all Sl-derived bands are below the 
Fermi level), so that Cu has valency 1.33, and the appro- 
priate ionic model is (Cu+''/3)3(S-2)(S-2). 



This means that the Cu d— band has 1/3 hole per Cu 
ion, 2.5 times fewer than in the high-Tc cuprates (optimal 
doping corresponds to 0.8-0.85 holes). This may be too 
far from half filling for strong correlation effects, but it is 
nevertheless suggestive of possible spin fluctuations. We 
will return to this point later. 



In order to understand the Cu d bands near the Fermi 
level, let us consider a simple tight binding model with 
two d orbitals, with m = ±2 (corresponding to combina- 
tions of the — and xy cubic harmonics, which be- 
long to the same representation in the hexagonal group). 
Since these orbitals are the ones spread most far in the 
plane, their hybridization with S is the strongest and they 
form the highest antibonding states near the F point, 
crossing the Fermi level. Integrating out the S Px,y or- 
bitals we arrive at the following model band structure: 
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where t are the Cu-S hopping amplitudes, and Ri are 
the three standard triangular lattice vectors, R^ — 0. 
Note that these bands are degenerate at F, unless spin- 
orbit coupling is taken into account. Near the top of the 
band the dispersion is isotropic, and away from it the 
Fermi surface develops a characteristic hexagonal rosette 
shape (Fig. g]). 

Let us now look at the calculated bands (Fig. [J). It 
is more instructive to concentrate on the righ hand side 
of Fig. [3J where the dispersion does not obscure the 
states degeneracy. We see, as predicted by the model, 
three sets of nearly parabolic bands, each four times de- 
generate at the point A=(0,0,7r/c). One of them is be- 
low the Fermi level and two above, forming the eight FS 
sheets we see in Fig. 21 The middle bands are predomi- 
nantly formed by the Cul and the lower (fully occupied) 
and upper ones by the Cu2, although there is substantial 
mixture of all three Cu orbitals. The average occupa- 
tion of Cu d orbitals, as described above, is 1/3 hole per 



Cu, too small to form a magnetic state, even in LDA-I-U 
with f7 ~ 8 eV (as verified by direct calculations). For- 
mation of an ordered magnetic state is additionally hin- 
dered by the fact that supexchange is this case is antifer- 
romagnetic, and frustrated, as it should be on a triangu- 
lar lattice. One may think that additional hole doping, 
achieved through Cu vacancies, broken S-S bonds or in- 
terstitial oxygen (note that this structures includes large 
pores, one per formula unit, in each Cu-S layer) should 
bring the c?— bands closer to half-occupancy and promote 
local magnetic moments. Note that in at least one exper- 
imental paper a Curie- Weiss behavior was reported, cor- 
responding to 0.28 /xs/Cu^, and in another a weak, but 
inconsistent with the Pauli law, temperature dependence 
was found^a, while other authors reported temperature- 
independent susceptibility. 

One can speculate that the unexpected high- 
temperature superconductivity observed by Raveau et 
al^ is a phenomenon of the same sort, namely that this 
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FIG. 4: (color online) Calculated Fermi surfaces in the hexag- 
onal (top) and orthorhombic (bottom) structure, viewed along 
the c-axis. Note reduced fcz dispersion in the bottom panel. 

superconductivity forms in a portion of a sample, the 
same portion where some previous researchers observed 
local magnetic moments. As discussed above, it is highly 
unlikely that a stoichiometric, defectless CuS sample 
would support either local moments or unconventional 
superconductivity. However, it is of interest to consider 
a hypothetical situation that would take place if such mo- 
ments were present. Indeed in that case one can write 
down the superexchange interaction between the nearest 
neighbors as antiferromagnetic Heisenberg exchange, in 
which case in the reciprocal space it will have the follow- 
ing functional form: 

J(q) =J^cos(k-R,). (2) 

i 

In Fig. [S]we show an example of a Fermi surface gen- 
erated for the model band structure (Eq. [1]), for the 
simplest case of tpd-rc — 0. The wave vectors correspond- 
ing to the peaks of the superexchange interaction ([2]) are 
shown by arrows. An interesting observation is that for 
this particular doping this superexchange interaction (or, 
better to say, spin fluctuations generated by this superex- 




FIG. 5: (color online) A model Fermi surface, calculated using 
Eq. [U overlapped with the wave vectors corresponding to su- 
perexchange on a triangular lattice. The signs show a possible 
/—wave pairing state, consistent with superexchnage-induced 
spin fluctuations. 

change) would be pairing for a triplet /—state shown in 
the same picture^. Indeed, the superexchange vectors 
always span the lobes of the order parameter with the 
same sign. Since in a triplet case spin-fluctuations gen- 
erate an attractive interaction, it will be pairing for the 
geometry shown in Fig. [51 Note that this is opposite to 
high-Tc cuprates, where the superexchange interaction 
spans parts with the opposite parts of the d— wave or- 
der parameter, but in a singlet channel this interaction 
is repulsive, and therefore pairing when the correspond- 
ing parts of the Fermi surface have opposite signs of the 
order parameter. 

While the model Fermi surface shown in Fig. [5] is 
roughly similar to that calculated in the stoichiometric 
CuS, the system at this doping is too far from the ordered 
magnetism to let us assume sizeable superexchange-like 
magnetic fluctuations. As mentioned, our attempts to 
stabilize an antiferromagnetic (more precisely, ferrimag- 
netic, since we only tried collinear magnetic patterns) 
using a triple unit cell failed, even in LDA-I-U. One may 
think of a hole doped system, where superexchange is op- 
erative and the inner Fermi surfaces (albeit not the outer 
ones) have geometry similar to that featured in Fig. [S] 

Of course, it may not be possible to stabilize a system 
at sufficient hole doping and retain the required crystal- 
lography. We prefer to think about the model discussed 
in the perevious paragraphs as inspired by the CuS cov- 
ellite, bit not necessary applicable to actual materials 
derived from this one. The reason we paid so much at- 
tention to it is that this is a simple generic model, describ- 
ing any triangular planar structure with transition metals 
and ligands in the same plane, as in the covellite, in case 
where correlations are sufficiently strong to bring about 
spin fluctuations controlled by superexchange. It is quite 
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exciting that, compared to the popular spin- fluctuation 
scenario of superconductivity in cuprates, to which it is 
conceptually so similar, this simple generic model results 
in a completely different superconducting state, triplet 



/, as opposed to singlet d. This finding may have im- 
plications far beyound this particular material and (yet 
unconfirmed) superconductivity in it. 
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